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ABSTRACT 

We Investigate the stability of electrostatic waves to 
thermodynamic and electric potential gradients. The major virtue 
of this analysis, other than Its overall generality. Is that 
thermodynamic gradients drive Instabilities even when the internal 
electric field vanishes. This result does not emerge from previous 
analyses because skewing of the distribution function was not 
Induvled In the dielectric. ^ 
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The role of lou-frequency electrostatic Instabilities In Inhibiting 
transport within laser produced plasmas has h^en the focus of substantial 
theoretical Irterest In the last few years. For example, Ion acoustic and 
electrostatic Ion cyclotlon Instabilities are known to be driven by non-thermal 
features of the electron distribution function associated with electric currents 
and/or heat fluxes, and can severely Inhibit transport within a plasma. The 
Importance of transport Inhibition In ^strophyslcal plasmas Is also being 
recognized - a nuoiber of articles addressing current sad heat flux limitations 
by electrostatic Instabilities during solar flares having also appeared In tho 
astrophyslcal llteraturo.^''^^ 

Since transport Inhibition In both laboratory and natural plasma systems 

appears to be of ixilversal Importance, a general linear analysis metlt? close 

attention. Despite numerous articles on this subject only Instahlllties 

directly dlrven by Inducti^ I electric fields have been considered. For exaniple, 

2 

Klndel and Kennel have Investigated the Ion acoustic and electrostatic Ion 

cyclotron Instabilities driven by resistive electric fields, B,reslstlve " 

In plasmas for which the net current Is non-zero. Another example Is the so- 

called "heat flw<" instabilities which /re driven solely by a cond>lnatlon of 

thermoelectric fields, ^thermo " II , and electron pressure electric 

fields, p “ ^11^^ , In plasmas for which the net current Is zero,^^"^^ 
me 

Under the zero current condition It Is noted In these articles that the 
Internal electric field, E||, adjusts such that E|| *■ - 

resulting In a non-thermal feature In the electron distribution function. This 

feature has been referred to as a "return current."^ There Is, of coursb, 

an additional non-thermal skewing of the electron distribution function directly 

related to the temperature and pressure gradients . However, unlike the 

E|| - field Induced "return current", the skewing of the distribution function 

by temperature and pressure gradients was not Included In the plasma dielectric. 
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Therefore, any resulting plasna instabilities are driven solely by Internal 
electric flelds->regardless of the source of these fields (eg. temperature 
and pressure gradients. Induction fields, etc.). Strictly speaking than, "heat 
flux" Instabilities, as addressed In the literature, r;re not categorically 
different from "current driven" Instabilities since the driving agent In both 
cases is an Internal electric field. 

In this note we Include, In the plasma dielectric, the skewing of the 

electron distribution function by temperature and pressure gradients and 

Investigate the resulting modifications to the stability threshold. It 

should be stressed chat retention of the nonthermal skewing allows us to 

study a class of Instabilities that are categorically different from those 

driven by internal electric fields, viz., temperature and pressure gradient 

Instabilities, which can occur even If the Internal electric field Is zero. 

Our analysis will also ^'elax the zero current constraint which a-orlorl 

defines the magnitude of the Internal electric field, B||. In. general, the 

internal electric field does not adjust such that j^| " 0. This Is especially 

* 

true during non-steady electrodynamic conditions when currents can be created 
by Induction. A determination of the self-consistent internal electric field 
Is highly model dependent, depending on a variety of sources of emf (eg. 
temperature and pressure gradients, suprathermal particle beams. Induction 
fields, etc.) and thus cannot ^e reasonably addressed In this note. There- 
fore, In our analysis, the Internal electric field will be explicitly treated 
as an undetermined parameter, along with the temperature and pressure gradients. 
In order to Isolate the different sources of Instability. The primary goal 
of this analysis Is to determine, as a function of the electron-ion temperature 
ratio, the threshold electric field, temperature, gradient, and pressure gradient 
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above which Ion-acoustic and/or electrostatic Ion cyclotron waves are 
destabilized. The major virtue of this analysis, other than Its overall 
generality, Is the realization that under certain conditional electron tem- 
perature and/or pressure gradients can drive plasma Instabilities even when 
the Internal electric ;!leld vanishes. This Important result does not emerge 
from previous analyses because the evaluated plasma dielectric did not 
include the gradient -induced skewing of the electron distribution function. 

The steady state electron distribution function, fe(v, Cg, cxi Cp), Is 
modeled by a Maxwellian plus small non-thermal components associated with the 

presence of an Internal electric field, cg ■ a temperature gradient, 

7 1 T VP 

«T " » a**** a pressure gradient, e_ ■ H 

e Sj 3 ^ ' 


feCv. eg* *T» *p^ ■ (rrvte^)"^ exp + E l^(v) cos 0)J (1) 

corresponding to the first two terms of a Legendre expansion in pitch angle, 9. 

In equation (I), Vg - (ZTg/mg)^, v^ - Vj, f + Vj.2, v„ « a - E, Tg, Pg, 

Ep Is the Drelcer field Eq *■ 4tt ne e^ InA/Xe and where: 

Dfey(v) - E (2) 

m e 

Is the analytic form In the weak anisotropy limit which follows from an 

14 

expansion of the distribution function In Sonlne polynomials (cf., Braginskll ). 

The coefficients Sq^ are obtained by a least squares fit to the results of 
15 ifi 

Cohen et al and Spltzer and Harm^” and are tabulated in Table I. Figure 1 
Illustrates the structure of fe for the two cases which isolate the non-thermal 
features associated with electric fields (le, eg ■ 0.15, ex * tp ■■ 0) and 
temperature gradients (le, e^ " 0.15, eg " ep " 0). Note that for each case 
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there la a region in velocity space for which ^^®/9v > 0, a necessary con- 
dition for lnstid>ility. This occurs for the higher velocity electrons. 
Furthermore, note that the current, J|| , and heat flux, Qp , are not 
necessarily zero, This follows from the first and third velocity moments 
of equation (1) which result in Onsager's relations, viz,, 

J„ - neVe <5.28 eg + H.26 tp + 3.70 ex) (3) 

3 o 

Q|l * - 7 oe me Ve'^(5y60 eg + 5.60 ep + 6,77 ex)* (4) 

Although the model electron distribution function used in this note does not 
Include: (1) Feedback from potentially excited plasma turbulence and 

(2) a runawiQT region in velocity space (l.e,, for e ”^), it is a 
reasonable choice for marginally stable systems provided « 1* 

The dispersion function for electrostatic waves in a « 1 magnetoplasma 
is: 


D<u), k) - 1 + E 

a-e,i ^ 


(5) 


where the electron and ion susceptibilities are given by: 


■ ;;2^T j ‘ [Eo(^ 

DG » A 


and 


Xi<ou,t) - 2 ^ S e-M-i I„(M.i) [1 + 5i 2(5^ - m?ci)l 

\kVi/ m"-« 


(7) 


with Xpg" 


4TTrae^ 


OU, 


'pi 


2 _ 4Tme2 
me 


, la - , ?ci 

k|| va 


^1 

k,|Vi’ 


2 gB 

■ %(hxVi/%) , , and where is the modified Bessel function, while 

Z(^) is the plasma dispersion function. In deriving these susceptibilities 
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aquation (1) hai baan uaad for tha alactron diatrlbution function and 
axp(-v^/v^^) for tha iona. Tha alactron auscaptlbllity 
coaf flcianti, An^ and in aquation (6) dapand upon tha alactron diatri- 
butlon througit tha a^ig^ in aquation (2)(c.f. Morriaon and lonaon^^) and aro 
tabulated in Table I. 

r 

Following tha uaual procedure of letting u) tu + if with ^ 0, tha 

dlaperaion function given by equationa (5)<-(7) reaulta in the following zero 
growth rate condition: 

‘ - s «« *00 - 5. [l * Z + 5,2] (8) 

Ihia condition, which follows from Im(D(u>,1()) * 0, correaponds to the zero 
growth rate relation between the magnitude of the electron distribution 
function's nonthermal features, characterized by e " E Cc^ Aoor, and marginally 
unstable waves of frequency u> and wavenumber 1^. In deriving equation (8) we 
have used « 1 thereby allowing us to define a generalized instability 
parameter, a “■•a convenient measure of the electron's tendency to drive electric 
field and pressure and temperature gradient instabilities. The absolute 
marginal stability condition is found by minimizing equation (8) for c with 
respect to and for fixed mi/me and Te/Ti. The resulting absolute 
marginal stability condition is given by: 
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where D, ■ 

Ti 

jen [2/ 


.3/2 n V * * 

i) rj . r„ - .-w 


Te 

®i L' 

ime/ \Tj 

L' j 


c* Is the minimum value o£ « above which th^ plasma is unstable to the electro- 
static waves and ~ 1.2 has been used^^. Note that equation (9) for 

** corresponds to an electrostatic ion cyclotron instability (which has the 
lower threshold for the temperature ratios shown) whereas equation (10) 
corresponds to an ion acoustic instability. These results are illustrated in 
Figure 2. 

The major emphasis in this Note has been to determine the form of a 

^ » 

general instability parameter, c "2.6 (c£'Hp) + 0.20 as a function cf 
Te where eg ■ , e.j ■ » and Cp ■ • These results reduce to those of 

Klndel and Kennel^^^ for *1 ■ «p “ 0. However, since Forslund^^ and Singer^^ did 
not Include gradient induced "skewing" of the electron distribution function in 
the dielectric, their results are somewhat different than ours. Specifically, 
for the zero net current condition, l.e.^^ cg + Cp " - .71 cx» they find thst 
*F.S. “ “ *T we find that a more accurate result is e* " - 1.65 ej. 

More importantly, however, our results indicate the possibility of 
instability even when the internal electric field is zero (i.e. instability 
can occur when eg ■ «p “ 0 at e* - 0.20 ex). This Interesting 
result only emerges when one includes the gradient induced skewing of the 
electrons In the dielectric and thus did not appear in previous analyses of 
this problem. 
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TABLE I, EUctcon dlitributlon function «nd tlactiron ■uicoptlblllty confflclents 
(c.f., aquations (1), (2), and (8)), 


m 


■»t 

^ibE"AbP 

^nT 



0 

0 

0 

2.6 

0.20 

-0.48 

1.0 

1 

-7.84 

4.75 

-2.6 

-0.20 

-7.6 

4.5 

2 

2.37 

-3.62 

8.4 

o 

a 

1 

2.3 

-3,5 

3 

-0.516 

0.626 

-2.5 

3.8 

-0.52 

0.62 

4 

5.72x10*2 

-7.67x10*2 

0.52 

-0.66 

5.6x10*2 

-8.8x10‘2 

5 

-3.17x10“3 

3.42x10*3 

-6.0x10*2 

7.8x10*2 

-3.2x10*3 

3.4x10*3 

6 

7.04x10*5 

-3.14x10*5 

3.2x10*3 

-3.4x10*3 

7.2x10-5 

-3.1x10-5 


-7.2x10“5 3.1x10"5 


7 







FIGURE GAFFIONS 

Figure 1. Electron diitrlbution function profiles es 4 function of the 
velocity component parellel to the magnetic field for several 
values of the perpendicular component. The upper curve assumea 
only an electric field is present, the lower that only a tempera- 
ture gradient exists. 

Figure 2. Instability threshold, as a function of electron to ion 
temperature ratio. 







